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ABSTRACT 

The thermal reprocessing hypothesis in AGN, where EUV/X-ray photons are repro- 
cessed by the accretion disc into optical/UV photons, predicts wavelength-dependent 
time delays between the optical continuum at different wavelengths. Recent photo- 
metric monitoring by Sergeev et al. has shown that the time-delay is observed in 14 
AGN, and generally seen to increase with increasing wavelength, as predicted in the 
reprocessing scenario. We fit the observed time delays and optical spectral energy dis- 
tribution using a disc reprocessing model. The model delivers estimates for the nuclear 
reddening, the product of black hole mass times accretion rate, and the distance to 
each object. However, the distances at face value give Hq = 44 ± 5 km s^^ Mpc^^ - 
a factor of 1.6 smaller than generally accepted. We discuss the implications of this on 
the reprocessing model. 

Key vifords: galaxies: active - galaxies: nuclei - galaxies: Seyfert 



1 INTRODUCTION 



Echo 



mappmR, 



reverberation 



mapping, 

l|Blandford fc McKed Il982l : iPetersonI Il993l ) uses light 
travel time to measure spatial separations within a dis- 
tant accretion flow. Ionizing photons produced close to 
the compact object irradiate the surrounding gas which 
reprocesses these into UV/optical continuum and emission 
lines. As the ionizing radiation varies erratically, so do 
the reprocessed components but with time delays due to 
light travel time within the system - the light travel time 
from source to reprocessing site to observer is longer than 
that on the direct path from source to observer. These 
observable delays provide indirect information on the size 
and structure of the surrounding accretion flows. The 
micro-arcsecond resolution that can be achieved by echo 
mapping is unreachable by any other technique (1 light day 
across at 100 Mpc is ~2 micro-arcsecond). 

It has long been known that the optical and UV con- 
tinuum in AGN are highly correlated and vary nearly si- 
multaneously. Initial attempts at measuring the time-delay 
(or lag) between the optical and UV continuum found i t 
to be less than a couple of days (e.g., IStirpe et al.lll993 ). 
However, higher time sampling by later spectroscopic mon- 
itoring campaigns have allowed for the measurement of 
wavelength-dependent continuum time delays in two objects 
(NGC 7469 and Ark 564). These two objects appear to fol- 
low the predicted relation for an irradiated accretion disc 
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l|Wanders et all Il997l : ICoUier et all Il998l . 1200 ll h with the 
lags being of the order of 1-2 days. Such small lags can only 
be explained in terms of AGN reprocessing as disc instabil- 
ities should produce time delays of the order of the viscous 
timescale (the timescale on which the local surface density 
changes) which is at least hundreds of d ays in AGN for th e 
region where the optical emission peaks (i Krolik et al.lll99ll ). 

Recent photometric monitoring by Sergeev et al.l (|2005f ) 
showed significant lags between optical lightcurves in 14 
AGN, with the lags between the lightcurves determined by 
standard cross-correlation techniques. Their measured lags 
range from tenths of a day to several days. The lag increases 
with wavelength, as is predicted by reprocessing in an ac- 
cretion disc, where the inner hotter regions see the ioniz- 
ing source before the outer cooler regions. They also find 
that the delay is systematically greater for greater absolute 
nuclear luminosity, following r x where b ~ 0.4 — 0.5, 
consistent with a disc reprocessing model. The lags are in- 
terpreted in terms of the light-travel time from an ioniz- 
ing source above the disc to the region in the disc where 
the ionizing radiation is reprocessed into optical continuum 
emission. Thus, the brighter the source, the greater the light- 
travel time to the region emitting the optical continuum. 

ICoUier et al.l (|l999D showed how measuring time delays 
between the optical/UV continuum in an AGN accretion 
disc and the fiux from the disc, the distance to an AGN 
can be measured, and hence AGN can be used as stan- 
dard candles to determine the current expansion rate of 
the Universe, Bubble's constant. Ho- From measurements of 
the wavelength-dependent time delays in NGC 7469, these 
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authors determine Ho = 42 ± 9 km Mpc^^ assuming 
an inclination of 45°. However, this is at odds with other 
determinations of Ho- Over the past decade, several meth- 
ods have been employed to determine Ho, with all methods 
now converginR to a v alue of around 72 km s~^ Mpc~^ (e.g. 
iFreedman et all 120011 ). The method using AGN accretion 
discs does not rely on the distance ladder calibration, and 
thus provides another independent check on these results. 
One distinct advantage of this method, compared to other 
methods, is that AGN are common in the Universe, and 
can be found over a wide range of redshifts. As they are 
some of the most energetic sources in the Universe, AGN 
can be seen at redshifts beyond the supernova horizon, and 
can thus potentially set tight constraints on ^Im and Qa - 
it is at redshifts greater than z ~ 1 where large differences 

are seen between different cosmologic al models. 

In this paper, we re-analyse the iSergeev et al.l (|2005| ) 
lightcurves, also in terms of thermal reprocessi ng in an AGN 
disc, i n an attempt to extend the method of ICollier et al.l 
j 19991 ) to use AGN accretion discs as a cosmological probe. 
In Section 2 we discuss the AGN lightcurves that we analyse 
and in Section 3 we detail the accretion disc model used. 
The results of our model fits to the data are presented in 
Section 4, and the implications of these results discussed in 
Section 5. 



2 AGN LIGHTCURVES 

We use the p u blishe d optical lightcurves for 14 AGN from 
ISergeev et al] l|2005l ) in the B, V, R, Rl and I bands. 
The fi lters used are non-standard filters l|Doroshenko et al] 
l2005al lb[). though their B, V, and R are close to Johnson fil- 
ters, the Rl is close to the Cousins I filter and the I filter 
is close to the Johnson I filter. The published lightcurves 
are given as relative fluxes with respective to comparison 
stars. However, the magnitudes of the comparison stars are 
only given in the standard B, V, R, I bands, rather than 
the non-standard bandpasses used. For this work, we re- 
quire the absolute fluxes of the AGN, rather than the rel- 
ative fluxes, and therefore require the magnitudes of the 
comparison stars in the filters used. To calculate these we 
use XCAL, Keith Home's interactive synthetic photometry 
program. XCAL computes predicted count rates from tar- 
get spectra in standard or non-standard passbands, allowing 
stellar spectra to be fit to the comparison star magnitudes 
in the standard Johnson- Cousins filters. Spectra of main 
seq uence stars from the B ruzual-Persson-Gunn-Stryker at- 
las (|Gunn fc Strvkedll983l ) were individually scaled to fit 
the observed magnitudes, and the star with lowest "^^^ 
adopted to model the spectrum of the comparison star. Af- 
ter finding the best-fit stellar spectrum, an intergral over the 
non-standard passband gives the magnitudes. In this way, we 
determined the fiuxes of th e comparison s tars i n the B, V, 
R, Rl and I filters used bv lSergeev et all (|2005l ) (see Table 
[T]and thereby convert the AGN lightcurves from differential 
magnitudes to absolute fluxes. 

All 14 AGN show significant variability over the ap- 
proximately 3 yea r monitoring period (see fig. 1 from 
ISergeev et al.|[2005h . The median temporal sa mpling for the 
lightc urves is 2 to 3 days (see table 3 from ISergeev et al.l 
I2OO5I . for sampling rates for each AGN). As an example of 



Table 1. Red shifts, z, and fluxes (in mjy) for the comparison 
stars used by ISereeev et al. for each of the 14 AGN in 

their B, V, R, Rl and I filters. 



Object 




z 


B 


V 


R 


Rl 


I 


NGC 4051 





002 


5.1 


7.9 


9.3 


9.7 


9.8 


NGC 4151 





003 


8.7 


26.6 


50.8 


66.4 


70.5 


NGC 3227 





004 


15.2 


28.5 


36.7 


39.5 


40.1 




n 
U 


uuy 


lU. 1 


10. D 


1 S 7 


1 Q A 




NGC 7469 





016 


26.5 


55.2 


77.2 


88.0 


90.2 


NGC 5548 





017 


6.0 


11.1 


14.4 


15.5 


15.7 


Mrk 6 





019 


3.5 


5.9 


7.5 


8.0 


8.1 


MCG-l-8-11-11 





020 


6.3 


19.2 


36.7 


48.0 


50.9 


Mrk 79 





022 


22.0 


37.6 


47.9 


51.0 


51.3 


Mrk 335 





026 


8.0 


13.6 


16.9 


18.1 


18.3 


Ark 120 





032 


4.5 


10.2 


15.0 


17.1 


17.5 


Mrk 509 





034 


25.0 


40.7 


49.2 


51.4 


51.8 


3C 390.3 





056 


3.8 


7.0 


9.0 


9.9 


10.0 


IE 0754.6+3928 





096 


5.8 


10.9 


14.0 


15.1 


15.3 



the strong correlations between the lightcurves in the dif- 
ferent wavebands, all 5 lightcurves for NGC 5548 and NGC 
7469, with the cross-correlation functions, are shown in Fig. 
[1] We have determined the lags between the lightcurves (rel- 
ative to the B-band lightcurve) for e ach AGN using the 
interpolation cross-correlation method l|Gaskell fc PetersonI 
ll987l : IWhite fc Petersonl F 



1994h. and find the d elays to be the 



Sergeev et al] HqoI) (see their Ta- 



same as determined by 
ble 4). The peak time delay is generally smaller than the 
centroid time delay, though in many cases they are consis- 
tent within the uncertainties. The maximum and minimum 
fiuxes of each of the lightcurves are also given in Table [2] as 
these are used in the model which we describe later. 



3 ACCRETION DISC MODEL 



ISergeev et al.l (|2005l ) successfully modelled the lag- 
luminosity relation in these sources with a disc reprocessing 
model, where interband time delays exist due to different 
light travel times from the ionizing source and the different 
continuum emitting regions. Here we take this further, 
by simultaneously fitting the interband lag distribution 
and spectral energy distribution (SED) in terms of a disc 
reprocessing model. In this model we fit the SED at the 
maximum and minimum points in the lightcurve, interpret- 
ing the change in brightness as a change in the luminosity 
of the ionizing source driving the reprocessing in the disc, 
and hence a change in disc temperature. Simultaneously we 
fit the lags with the same model. 

We assume that there is a highly variable source of ion- 
izing radiation near the disc axis that illuminates the disc, 
driving optical/UV variability and leading to wavelength- 
dependent time delays. The hot inner regions which emit 
mainly UV photons 'see' the driving ionizing photons be- 
fore the cool outer regions which emit mainly optical pho- 
tons. Thermal radiation from a disc annulus at tempera- 
ture T{R) emerges with a range of wavelengths, centred at 
A ~ hc/kT{R). Roughly speaking, each wavelength picks out 
a different temperature zone and the time delay t = R/c 
measures the corresponding radius. Thus, shorter wave- 
lengths sense disc annuli at higher temperatures. 
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NGC5548 NGC7469 




2400 2600 2800 3000 -20 20 2600 2800 3000 -20 20 

JD - 2450000 T (days) JD - 2450000 t (days) 



Figure 1. Lightcurves for NGC 5548 and NGC 7469 in all 5 bands. Corresponding cross-correlation functions (CCF) for the lightcurves 
with respect to the B-band are also shown. In the B-band the auto-correlation function is plotted. Solid line indicates a time-delay, r = 
0.0 d, and the dotted line shows the centroid of the CCF. 



In this model there is heating due to irradiation as well 
as viscous heating in the disc. For viscous heating alone, 
when the radius is much greater than the inner most stable 
orbit, the disc temperature profile should follow 



T(i?) = 



/ 3GMM \ 



1/4 



(1) 



where M is the mass of the black hole, and M the mass 
accretion rate. In the case of irradiation of the disc from 
a height above the disc, the temperature profile is a 
combination of the temperature due to irradiation and the 
temperature due to viscous heating, 

1/4 



T{R) 



f sgmm \ 



{1-A)L, 



(2) 



where Lx is the luminosity of the irradiating source, A is 
the disc albedo, i?* is the distance from the illuminating 
source to the surface element of the disc, and dx is the angle 
between the disc surface normal and a unit vector from the 
surface element to the variable central source. This leads to 
T{R) oc (LxHxf'*R-^'*, for R > Hx and provided th e 
disc thickness is H ^ Hx (e.g. iFrank. King fc Rain"3l2002l ). 
Thus, in both cases, we expect T{R) oc J?"^''*. As the time 
delay goes as, t = R/ c and the wavelength goes like A oc , 
then this predicts that the wavelength-dependent time-delay 
profile follows a r oc A^''^ relationship. 

With the temperature profile determined from the 
wavelength-dependent time delays, the predicted spectrum 
can be determined by summing up blackbodies over disc an- 
nuli, and comparing this with the observed spec trum leads to 
a det ermination of the distance to the source (|Collier et al.l 
119991 ). In this way the predicted spectrum is given by 



2hc 



2nR d_Rcos i 



\Z ghc/AfcT . 



D2 



(3) 



with the Planck function, AO, the solid angle subtended 
by the disc annuli, Rin and Rout the inner and outer disc 



radii, i th e inclination a ngle of the disc, and D the distance. 
Following IColheil I 19981 ) ■ we assume a general temperature 
profile for the disc of T = Toy'*' with y = R/Ro {To is 
the temperature at radius Ro), and so — hc/XkTo such that 
s = Soy'' and Ay = y As/bs. Over a range of wavelengths 
that are well within those corresponding to the inner and 
outer disc radii, substituting these expressions into Eq. [3] 
leads to 



47ri?o cos i 



D2 



(/ic)^-'^/'"A(^/'')-^(fcro)^/'':^ 



where 



/„(&) = 



5{n/6)-l 



(4) 



(5) 



7„(fo) is a numerical function, and we calculate J2{3/4) = 
1.93 and /3(3/4) = 6.49 (6 = 3/4 corresponds to the stan- 
dard thin disc assumption). 

AGN spectra, however, rarely match the classical thin- 
disc spectrum oc i/^^^. This difference is likely due to 
contributions from stars and dust. Nevertheless, the disc 
spectrum can be is olated by taking difference spectra, and 
ICollier et al.l ()l999[ ) found that the variable component of 
AGN light in NGC 7469 does have A/^ oc u^^^. For this rea- 
son we choose to model the difference spectra of the AGN 
using the method described below. 

3.1 Disc Transfer Function 

We model the AGN wavelength-dependent time-delay distri- 
bution and spectrum using a disc transfer function 'I'y(r, A) 
to describe the response of the disc to cha nges at a give n 
wavelength and time delay, as outlined by IColheil l|l998h . 
The disc transfer function is defined as 



dB, dT 
dT dLx 



R 

S(t (1-1- sin i cos 9)) dfi 

c 



(6) 
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Table 2. Maximum and minimum fluxes (mjy) from the AGN 
lightcurvcs. Continued on next page. 



Table 2. continued. 



Object 


Filter 


Maximum flux 


Minimum flux 






(mJy) 


(mJy) 


NGC 4051 


B 


12.89 ±0.12 


9.92 ± 0.05 




V 


21.46 ±0.16 


18.22 ±0.08 




R 


34.23 ±0.19 


29.71 ±0.16 




Rl 


41.31 ±0.27 


36.55 ±0.20 




I 


43.42 ± 0.30 


38.94 ±0.24 


NOG 4151 


B 


57.53 ±0.54 


29.93 ±0.29 




V 


74.97 ±0.35 


49.07 ±0.27 




R 


117.93 ±0.49 


84.37 ± 0.34 




Rl 


125.98 ±0.54 


91.60 ±0.38 




I 


133.10 ±0.56 


97.44 ±0.43 


NGG 3227 


B 


11.41 ±0.07 


8.69 ±0.05 




V 


22.84 ±0.13 


19.53 ±0.10 




R 


41.02 ±0.20 


36.95 ±0.18 




Rl 


51.03 ±0.26 


46.72 ±0.23 




I 


54.54 ±0.27 


50.11 ±0.52 


NGG 3516 


B 


15.85 ±0.10 


12.03 ±0.06 




V 


31.00 ±0.13 


27.33 ±0.13 




R 


52.15 ±0.22 


46.25 ±0.20 




Rl 


62.06 ± 0.39 


56.58 ±0.46 




I 


65.25 ±0.36 


59.74 ±0.30 


NGG 7469 


B 


12.28 ±0.06 


9.72 ± 0.04 




V 


19.53 ±0.09 


17.16 ±0.07 




R 


33.95 ±0.14 


30.73 ±0.13 




Rl 


41.56 ±0.17 


38.33 ±0.16 




I 


44.99 ±0.18 


41.76 ±0.17 


NGG 5548 


B 


6.94 ±0.05 


4.15 ±0.02 




V 


12.20 ±0.19 


9.33 ±0.04 




R 


20.17 ±0.11 


15.73 ±0.09 




Rl 


23.24 ±0.21 


19.09 ±0.13 




I 


24.71 ±0.10 


20.71 ±0.15 


Mrk 6 


B 


4.50 ±0.03 


2.99 ± 0.02 




V 


9.42 ± 0.06 


7.25 ± 0.04 




R 


17.05 ± 0.09 


13.45 ±0.06 




Rl 


19.07 ±0.11 


15.50 ±0.08 




I 


20.26 ±0.10 


16.44 ±0.09 


MGG+8-11-11 


B 


6.83 ±0.04 


3.61 ±0.03 




V 


10.60 ±0.05 


7.08 ± 0.05 




R 


19.61 ±0.08 


14.98 ±0.07 




Rl 


22.51 ±0.09 


17.53 ±0.08 




I 


24.44 ±0.10 


19.16 ±0.09 


Mrk 79 


B 


6.74 ±0.05 


4.02 ± 0.02 




V 


9.65 ± 0.08 


6.87 ±0.04 




R 


15.12 ±0.09 


11.50 ±0.05 




Rl 


16.22 ±0.12 


12.67 ±0.06 




I 


16.96 ±0.13 


13.32 ±0.08 


Mrk 335 


B 


7.61 ±0.06 


4.45 ± 0.02 




V 


8.90 ± 0.06 


6.04 ± 0.04 




R 


12.67 ±0.07 


9.00 ± 0.04 




Rl 


12.09 ±0.07 


8.74 ±0.04 




I 


12.58 ±0.08 


9.10 ±0.05 






-Lz.uo in U.-LO 


7 n'^9 -t- n 1 n 




V 


16.37 ±0.13 


10.60 ±0.08 




R 


23.55 ±0.13 


16.00 ±0.09 




Rl 


24.02 ±0.15 


17.04 ±0.09 




I 


25.02 ±0.16 


18.06 ±0.10 


Mrk 509 


B 


15.79 ±0.11 


12.11 ±0.09 




V 


18.07 ±0.11 


15.02 ±0.11 




R 


25.40 ±0.12 


22.24 ±0.11 




Rl 


23.57 ±0.13 


20.76 ±0.12 




I 


24.47 ±0.14 


21.71 ±0.12 



Object 


Filter 


Maximum flux 


Minimum flux 






(mJy) 


(mJy) 


3G 390.3 


B 


2.86 ± 0.02 


1.05 ± 0.01 




V 


4.12 ± 0.02 


2.19 ± 0.01 




R 


6.30 ± 0.03 


3.74 ±0.02 




Rl 


6.77 ± 0.04 


4.25 ±0.03 




I 


7.13 ±0.05 


4.53 ±0.03 


IE 0754.6+3928 


B 


5.12 ± 0.03 


4.45 ± 0.03 




V 


6.87 ± 0.04 


6.33 ±0.06 




R 


8.33 ±0.05 


7.62 ± 0.04 




Rl 


9.40 ±0.06 


8.58 ± 0.05 




I 


8.81 ± 0.05 


8.11 ±0.06 



0.8 



0.6 



0.4 

s 

-> 0.2 




10 



Figure 2. Irradiated accretion disc transfer functions for 4400A 
with disc inclinations of 0° (solid) and 45° (dashed) and 8OOOA 
with disc inclinations of 0° (dotted) and 45° (dash-dot). Tb = 
15000 K and Tp = 12000 K. Gircles mark the centroid delay for 
comparison with the peaks. 



-Rout 



dB^ dT R AR AO cos i 



dT dL^ 



D2 



&{t (1 + sin i cos 0)) 



(7) 



where is the driving ionizing luminosity. The 5-function 
ensures that only radii corresponding to the specific time 
delay contribute to the transfer function. As we are fitting 
the maximum and minimum SED of the AGN with this 
model, we treat the disc as having temperature Tb at Rq 
during the bright state and Tp at Ro during the faint state. 
Hence, 



dB^ dT 
dT dL^ 



= B,(T2)-B,(Ti 



(8) 



and T2 = Ts (-R/-Ro)"^ Ti = Tf{R/Ro)~'', and we take 
_Ro = 1 light-day. Examples of the transfer function are given 
in Fig. H 
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3.2 Time delays 

The centroid of the transfer function is equivalent to 
the centroid of the CCF which is used as a measure of 
the luminosity -weighted radius of the reprocessing region 



Koratkar fc Gaskel^ 1991.1 . however, see [Robinson & Perez 
1990l ) and lWelshI (| 19991) for limitations. For our blackbody 



disk model, the time delay centroid of the transfer function 



(r) 



J rVE'^(r, A) d? 



/ *4r,A) dr ' 
which becomes 

_ /rr'[g4r2) - B,iT,)]{Ryc) dR 
/iJr' [^-(^2) - B4Ti)]R dR 
This reduces to 



(9) 



(10) 



(11) 



_ Ro / XkTB y^" /3(b) 1-e'/' 
c \ he J h{h) 1-e 

where e = {Tf /Tb)'^^*' ■ Rearranging the above equation for 
fcTs and substituting into Eq. |3J the bright state disc spec- 
trum is given in terms of the measured lag; 

h{hf 



B 47rCOSi 3\-3/ \2 J-2 

= n-, he X (r) 



D2 



bhib)' 



( V 



(12) 



Similarly, it can be shown that the faint state spectrum is 
just given by — ef^ and the difference spectrum between 
the bright and faint state is given by A/j/ = ^ ~ (1 — 
e)/if • As here we will always use b — 3/4 (as is appropriate 
for a classical thin accretion disc), we evaluate Eg. 1111 and 
[IHfor & = 3/4 below: 



0.62 d 



_JI_ 
40.1 Jy 



(io^a) 



V io^k) 



4/3 



-3/2 



(13) 



lo^A 



7— y 

\^Mpc ) 



l-£ 
- £3/2 



(14) 



Therefore, if the bright and faint state fluxes and interband 
time delays are measured from monitoring observations of 
an AGN, then the distance to the AGN can be determined: 



D 



i.3Mpc 



A \-^^^ f / cosi 



lO^A 



Jy 



-1/2 



(1-53/2) 



(15) 



Thus from the reverberation results we obtain a distance 
estimate that is independent of reliance on the Hubble law, 
and so can be used to test cosmological predictions for the 
luminosity distance versus redshift relationship (see Section 
4.1 and 5). We also note that in the above equation all fac- 
tors of (1 + 2) to change from observed to emitted quantities 
cancel out. 



3.3 Applying this model to the Sergeev et al. data 

With interband time dela ys and long- t erm monitoring 
lightcurves for 14 AGN, the lSergeev et al.l l|2005l ') data offer 



an opportunity to test this reprocessing model and to deter- 
mine distances to a sample of 14 AGN. We fit the difference 
spectra (between the maximum and minimum points on the 
lightcurves) with the above accretion disc model. An ad- 
ditional component to the model not discussed above is ex- 
tinction due to interstellar reddening within our own Galaxy 
and also the within the AGN - this is required to correct the 
difference spectra to show a u^^^ slope. Therefore, as part 
of the model, we include E{B ~ V) as a, free parameter (we 
discuss this assumption below). The other free parameters 
in the fit to the data are the bright state temperature, Tb, 
the faint state temperature, Tf, the distance to the AGN, 
D, and the host galaxy fluxes in the V, R, Rl and I bands. 
For a lower limit on the distance (and hence an upper limit 
on Ho), we set the host galaxy contribution in the B band to 
zero. Whilst this is clearly an incorrect assumption, without 
it, a fit would be degenerate since a fainter accretion disc can 
be offset by a brighter galaxy. Independent measurements 
of the host galaxy in the bands, would constrain these nui- 
sance parameters to reduce or remove this degeneracy. In the 
model we assume that the accretion disc inclination angle, 
i = 45°, as would be expected for type I AGN, where the 
broa d emission lines a re visible, from unified AGN theories 
(e.g., lAntonucci|[l993 ). However, the affect on the distance 
determined is small as Z? ~ (cosi)'^/^. 

3.4 Extinction within the central regions of AGN 

Two of the assumptions made when applying this model are 
that the discrepancy in slope between the difference spectra 
and a u^''^ slope is due to extinction from our own Galaxy 
and within the central regions of the AGN, and also that 
the variable component of the AGN spectra does not change 
shape between the bright and the faint state. To test these 
assumptions, we use the flux varia t ion gr adient method of 
IWinkler et all l|l992l) and IWinklerl (119971 ). In this method, 
when comparing the flux through one filter with the fiux 
through another filter at corresponding times, a very strong 
linear relationship is seen. We show a couple of examples of 
'fiux-flux' diagrams in Fig. [3] obtained with the Sergeev et 
al. data. We choose to show flux-flux diagrams for the two 
most variable AGN in the sample (3C 390.3 and NGC 4151). 
Firstly, as the relationship is so close to linear, it demon- 
strates that the variable component has a virtually constant 
flux distribution, thereby confirming our assumption that 
this is the case. If the spectrum did, in fact, change shape, 
the relationship in the fiux-fiux diagram would curve rather 
than be linear. A lag between the continuum bands of a few 
days only introduces a small scatter about the straight line 
in the fiux-fiux diagrams, and our diagrams are consistent 
with this. 

The slopes of the best-fit lines give the ratios of the vari- 
able component of the fiuxes and hence the SED of the inner 
nucleus, as viewed from outside that galaxy. Of course, the 
intrinsic SED of this region will differ due to extinction from 
dust in the outer parts of the nucleus. Assuming some stan- 
dard colours for the nuclear SED, the extinction can then 
be calculated by comparing the standard colours with the 
colours determined from the fiux-fiux diagrams. Initially, we 
correct for t he extinction in our Galaxy, E{B — V)g, using 
values from [Schlegel. Finkbeiner fc Davii (| 19981 ). Next, we 
correct for extinction within the nucleus of the AGN. As- 
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NGC4151 3C390.3 




V flux (mJy) V flux (mJy) 

Figure 3. Flux-flux diagrams for NGC 4151 (left) and 3C 390.3 (right). B (squares), R (filled circles), Rl (open circles) and I (triangles) 
flux is plotted against V flux at corresponding times. The best straight line fits are also shown. NGC 4151 and 3C 390.3 are shown as 
these objects display the largest variability. 



suming that the intrinsic SED for the nuclear region of an 
AGN follows fiy oc f^^^ in the optical, as expected for a 
classical thin accretion disc, we adjust the nuclear extinc- 
tion, E{B — V)n, so that the reddened model matches the 
observations. 

This, of course, assumes that the difference between the 
spectra of the variable components and a u^^^ spectrum is 
due to extinction. As an independent check of the values we 
get from the flux variation gradie nt method, we also us e the 
Balmer decrement method (e.g. [Reynolds et al.lll997^ . As 
extinction by cosmic dust is highly wavelength dependent 
it can change observed line flux ratios significantly away 
from intrinsic values. In this method, the observed Ha to 
H/3 broad line ratio is compa red with an assumed int rinsic 
ratio, which we take to be 3.1 (jGaskell fc Ferlandll 19841 ) with 
the difference between the two determining the extinction re- 
quired. We use values in the literature for (Ha +[NII])/H/3 
and H/3 /[OIII]A5007, wh ic h are listed in a ra n ge of p revious 
studies: lAndersonI lll970l): iBoksenberg et al] (Il975l): [Coheiil 
il983l ):lde Griip et al' ('l992^:'Koski' ('l978'):'M orris fc WardI 



Ostcrbrock & Koski (1976); Ostcrbrocl 



1977 ) 



PhiU ips (197S): iRafanelli fc Schulzl (|199ll ): IStephensl l|l989D : 
Winkler, (1992h . From these measured ratios we determine 



the Ha to H/3 broad line ratio assuming that the broad 
and narrow components are similarly reddened. We do 
this by assuming that H/3(narrow)/[OIII]A5007 = 0.1 and 
Ha(narrow)/[NII] — 1, which are the approximate values 
where th ese ratios cluster for Seyfert 2 galaxies in the dia- 
grams of lVeilleux fc Qsterbrockl (Il987l ). 

Using these two methods the nuclear extinctions have 
been calculated, and the results are presented in Table |31 
A comparison of the reddening determined using these two 
methods is shown in Fig. U Whilst the nuclear extinc- 
tions determined by these two methods do not agree exactly 
(largely due to observational difficulties in obtaining accu- 
rate spectrophotometric data by the authors of the various 
studies listed earlier), those AGN with higher nuclear ex- 
tinction from one method also have high nuclear extinction 



Table 3. Nuclear extinction, E{B — V)n, determined by the flux- 
flux and Balmer decrement methods, as well as the extinction due 
to our own Galaxy, E{B — V)g. The total extinction is just the 
sum of both components. 



Object 


E{B - V)g 


E{B - V)n 


E{B - V)„ 






(flux-flux) 


(Balmer decrement) 


NGC 4051 


0.013 


0.17 


0.21 


NGC 4151 


0.028 


0.15 


0.04 


NGC 3227 


0.023 


0.26 


0.31 


NGC 3516 


0.042 


0.16 


0.15 


NGC 7469 


0.069 


0.04 


0.09 


NGC 5548 


0.020 


0.16 


0.28 


Mrk 6 


0.136 


0.28 


0.48 


MCG-(-8-ll-ll 


0.217 


0.13 


0.34 


Mrk 79 


0.071 


0.09 


0.20 


Mrk 335 


0.035 


0.07 


0.00 


Ark 120 


0.128 


0.00 


0.04 


Mrk 509 


0.057 


0.00 


0.11 


3C 390.3 


0.071 


0.14 


0.26 


IE 0754.6-f3928 


0.066 


0.00 


0.04 



from the other, and similarly for low values (see Fig.|4]|. This 
confirms that there is significant extinction in the central re- 
gions of these AGN that needs to be accounted for in our 
models. 



4 RESULTS 

We use the model described in the previous section to fit the 
data for all 14 AGN from the Sergeev et al. sample. Specif- 
ically, we fit the model to the time delays, and the maxi- 
mum and minimum SED by minimising using a downhill 
simplex (the 'am oeba' from Numerical Recipes in Fortran, 
iPress et al.lll992l) to determine the best-fitting non-linear 
parameters, whilst optimising the linear parameters using 
optimal scaling. For the time delays we adopt the peak of 
the CCF rather than the centroid. The peak time delays 
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Figure 4. Comparison of E(I} — V)„ determined by the flux- 
fiux and Balmer decrement methods. The dashed line indicates a 
one-to-one relationship. 
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Figur e 5. Solid line: the adopted extinction law of lNandv et al.l 
lIlOTSl l and ISeatonI l|l979l ') evaluated at E( B - V) = 1.0 for R = 
3.2. Dashed line: the AGN extinction law of|Gaskelletal. || |2004) 
which has R = 5.15. We have extended the G askell et al.l l|2004l ') 
law to longer wavelengths (determined up to A = 6250A). We 
use the lNan dv et al. (1975) and Scaton (197S|) law to correct for 
the known reddenin g in our Galaxy, and the extinction law of 
I Gaskell et al. I | |2004) to correct for the reddening intrinsic to the 
AGN. 



are generally smaller than the centroid delays and there- 
fore gives a lower measurement of the distance. Firstly we 
correct for the known galactic extinction in the direction 
of each AGN using the c ombined galac tic extinction law of 
iNandv et al] ()l975l ) and ISeatonI (j 19791 ). Further reddening 
due to extinction intrinsic to t he AGN i s inclu ded using the 
AGN reddening curve of .Gaskell et al.l l|2004l ) (see Fig. [S]). 
While these reddening laws both have very similar slopes 
in the optical, they differ greatly in the UV (where we are 
not concerned about here), and importantly, have different 
R values, where R = Av/E{B - V). 



Figure 6. error parabola for the distance to NGC 7469. Solid 

lino is the best-fit parabola to the values determined at a range 

of distances either side of the best-fit distance. The dotted line 

indicates the location of the best-fit distance. Dashed horizontal 

lines indicate . and . - 

l-CT errors in the distance, i.e., 
2 



- 1. Dashed vertical lines indicate 
the distances that correspond to 



X 



+ 1. 



To determine the uncertainties in each of the fit param- 
eters, we fix the parameter of interest at a value offset from 
the best-fit value. The other parameters are then optimised 
and the value determined. We do this for a range of val- 
ues either side of the best-fit value, with the range chosen 
so that it covers Ax = 1-0. A parabola is then fit to the 
points, and forced to be a minimum (and hence have zero 
gradient) at the best-fit value. From the parabola 1-cr errors 
are determined from where Ax^ = 1-0. An example of this 
is shown in Fig. (6) 

We present the best-fitting models for each of the AGN 
in Fig. [3 and the best-fitting parameters in Table U] Note 
that the extinction values are close to those determined by 
the two other methods (see Fig. |8] for a comparison), fn sev- 
eral cases there are deviations of the difference spectrum 
from the best-fit model resulting in high values of reduced 
X^- This is likely to be due to contributions from broad emis- 
sion lines. For instance. Ha lies in the R-band, and various 
Balmer lines and Fe If features lie in the B-band (see Fig. 
^ . The R-band fiux does sometimes appear higher than the 
model fit, particularly in cases where the x^ is high, as would 
be expected due to contributions by Ha. 

To check whether the determined galaxy spectra are 
realistic, we plot the galaxy spectra dereddened for galac- 
tic interstellar extinction, all normalised to the I band fiux 
in Fig. 1101 We do not deredden for the nuclear extinction 
that we have determined, as this reddening is expected to 
be due to dust in the nuclear region and not the galaxy as 
a whole. Additionally, we plot the galaxy spectrum for a 
Sab galaxy (as expecte d for a Seyfert ) using the intrinsic 
galaxy colours of .Fukugita et al.l (Il995l ). Whilst the major- 
ity of galaxy fiuxes match the model galaxy in the V and Rl 
bands, all of the galaxies have a higher flux than the model 
in the R band. We suggest that this is likely due to strong 
Ha line emission from the AGN that falls within this band 
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Figure 8. A comparison of the total E(B-V) values determined 
via the reprocessing model and the flux-flux and Balmer decre- 
ment methods. Circles indicate points determined by the flux-flux 
method, and crosses indicate those determined by the Balmer 
decrement method. Dashed line indicates a one-to-one relation- 
ship. 



(see Fig. I^}, making the R band galaxy flux higher than for 
a typical Sab galaxy. To check this we load the spectrum 
of 3C 390.3 shown in Fig. |9]into the synthetic photometry 
package XCAL, and calculate the predicted flux in the R 
band. We then subtract the Ha line from this spectrum and 
calculate the predicted flux in the R band once again. We 
find a factor of 1.3 difference between the 2 spectra. Com- 
paring the 3C 390.3 R band galaxy flux that we determine 
with that expected for an Sab galaxy we flnd a factor of 1.5 
discrepancy. Given that the Sab galaxy model is not going 
to perfectly match the actual galaxy spectrum of 3C 390.3, 
it seems reasonable that the observed difference could be 
due to Ha. 

There is one galaxy spectrum, IE 0754.6-1-3928, that is 
significantly different from all the others. In particular, the 
Rl flux is much higher than the other galaxies. This object 
has the highest redshift of the sample, with z = 0.096. The 
Ha line will have therefore shifted from an emitted wave- 
length of 6563 A to 7193 A, and into the Rl fiher. Note that 
the next highest redshift object in the sample is 3C 390.3, 
at z = 0.056, not high enough to shift the line significantly 
into the Rl band (see Fig. |9j . 



4.1 The Hubble Constant 

Having measured distances to the 14 AGN using the con- 
tinuum reverberation method detailed above, we are now 
in a position to determine Hubble's constant. Ho- Assum- 
ing a cosmological constant form for Qa, in the Friedmann- 
Robertson- Walker metric, the luminosity distance is deter- 
mined by: 



c 



-(i+z)ififcr^/'5fc|infcj 




5000 6000 7000 8000 

Wavelength (A) 

Figure 9. Optical spectrum of 3C 390.3, with the filter pass- 
bands shown (dashed line). The filters (from left to right) are 
the Crimean B, V, R, Rl and I. Note that Hq sits in the R- 
band. 3C 390.3 is the object with the second highest redshift in 
our sample (z = 0.056). Hence, in all other objects (expect IE 
0754.6+3928, which has z = 0.096), the Ha line will be closer to 
its emitted wavelength of 6563 A. Data from AGN Watch website 
( http: / /www.astronomy.ohio-state.edu'-^agnwatch/ [ I. 




5000 



6000 7000 
Wavelength (A) 



8000 



Figure 10. Dereddened galaxy spectra normalised to the I band 
(circles). Note the B-band galaxy flux is forced to be zero in the 
fit. The crosses show the galax y spectrum for a n Sab galaxy using 
the intrinsic galaxy colours of .Fukugita et al ] 1(1995 '). 



1+2 



dx 



\l/2 



(16) 



1/2 



where x = 1 + z, flk — 1 — — SIa and Sk is sinh for Qk > 
and sin for ilk < 0. Fo r Qk = 0, Eg. I16lrev erts to cHq^{1 + 
z) times the integral (|Carroll et al.l 19921 ). We assume the 
current standard cosmology with Qm = 0.3 and Qa ~ 0.7 
when fitting this model. As there is quite a large scatter 
in the distances (see Fig. Illjl . we also include a fractional 
systematic error, /o, as a parameter in the fit to account for 
this intrinsic scatter. The uncertainty on a single distance 
measurement, D, is then given by a-f{D) — a'j^ + /qD^, 
where is the measurement uncertainty on the distance 
from the fit. One cannot minimise when including such a 
term in the fit as continuing to increase /o will continue to 
decrease with no end. Instead we maximise the likelihood, 
which is equivalent to minimising — 21nL, which includes a 
penalty for expanding the error bars: 
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Figure 7. Reprocessing model fits for all 14 AGN. Each double panel (described next) shows the fit for one AGN. Top: Maximum 
(circles) and minimum (triangles) spectra, with best-fitting model. Difference spectrum (crosses), with reddened accretion disc model. 

De-reddened disk model is given by the dashed line. Squares indicate best-fit galaxy fluxes. Bottom: Time-delay distribution with best 
fit model. Time delays relative to B-band, and determined from the peak of the CCF. Figures continued on the next page. 
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Figure 7. continued. 

Table 4. Best-fitting parameters from reprocessing model. MM is calculated using the faint state temperature. This gives a maximum 
value for this product, as the temperature in the faint state will not entirely be due to viscous heating and there will be some contribution 
from irradiation. 



Object z D cz/D E{B-V)„ Tp Tb MM xl 

(Mpc) (km s-i Mpc-i) (10^ K) (10=^ K) (lO'^ yr-i) 



NGC4051 


0.002 


19 ± 17 


37 ± 33 


0.24 


± 


0.02 


4.1 ± 2.7 


4.5 ± 3.0 


0.03 


0.8 


NGC4151 


0.003 


19 ± 2 


52 ± 5 


0.17 


± 


0.01 


5.4 ± 0.5 


6.8 ± 0.6 


0.09 


8.1 


NGC3227 


0.004 


31 ± 10 


37 ± 12 


0.24 


± 


0.02 


5.7 ± 1.3 


6.3 ± 1.5 


0.1 


0.2 


NGC3516 


0.009 


63 ± 27 


42 ± 18 


0.20 


± 


0.02 


10.3 ± 3.3 


11.4 ± 3.6 


1.1 


3.5 


NGC7469 


0.016 


77 ± 22 


64 ± 18 


0.11 


± 


0.02 


9.6 ± 2.0 


10.4 ± 2.1 


0.8 


2.3 


NGC5548 


0.017 


341 ± 62 


15 ± 3 


0.22 


± 


0.01 


24.6 ± 3.3 


29.9 ± 4.0 


36 


5.9 


Mrk6 


0.019 


95 ± 30 


59 ± 19 


0.30 


± 


0.01 


11.8 ± 2.8 


13.8 ± 3.3 


1.9 


5.7 


MCG+S-ll-ll 


0.020 


248 ± 49 


25 ± 5 


0.08 


± 


0.01 


18.4 ± 2.7 


23.3 ± 3.4 


12 


3.9 


Mrk79 


0.022 


180 ± 43 


37 ± 9 


0.13 


± 


0.01 


13.5 ± 2.4 


16.4 ± 2.9 


3.3 


3.0 


Mrk335 


0.026 


219 ± 53 


35 ± 9 


0.10 


± 


0.01 


14.4 ± 2.6 


17.5 ± 3.2 


4.2 


7.3 


Arkl20 


0.032 


190 ± 51 


51 ± 14 


0.08 


± 


0.01 


16.8 ± 3.3 


20.9 ± 4.1 


7.9 


5.7 


Mrk509 


0.034 


376 ± 89 


27 ± 6 


0.00 


± 


0.00 


26.5 ± 4.7 


29.0 ± 5.1 


48 


2.6 


3C390.3 


0.056 


400 ± 205 


42 ± 22 


0.16 


± 


0.01 


15.7 ± 6.1 


22.8 ± 8.8 


6.0 


7.1 


IE0754.6+3928 


0.096 


1289 ± 130 


22 ± 2 


0.07 


± 


0.03 


54.3 ± 9.4 


57.0 ± 4.6 


856 


1.0 



N 

-21nL = x^ + 2^1n(cri)+const , (17) 

i=l 



where there are data points, and L ~ 

exp(— X /2)/{Y[ - V^TTai). The resuhing fit is shown in 
Fig.[TTl We find,' Ho = 44 ±5 km s"^ Mpc"^ and fo = 0.30, 
with a reduced = 0.91. The implications of this small 
value of Ho are discussed in the next section. 
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To check the possible range in distance if the B-band 
galaxy fiux is allowed to be non-zero (which it clearly will 
be), we also fit the disc reprocessing model when allowing 
the B-band galaxy fiux to be the maximum possible value, 
i.e., equal to the minimum flux from the lightcurve. Such 
a fit gives the maximum distance (and hence the minimum 
Ho) possible. Fitting to these distances gives Ho = 21±2 km 
s~^ Mpc~^. The distances and best-fitting model are shown 
in Fig. 111! as open circles and the dotted line, respectively. 
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Figure 11. Left: Hubble diagram for 14 AGN from the Sergeev et al. sample. Distance modulus is plotted vs redshift. The solid line is 
the best-fitting model (fit to the filled circles), where as the dashed line indicates Hq = 72 km s~^ Mpc~^. The filled circles are for the 
distances determined when the B-band galaxy flux is = 0, and the open circles are for when the B-band galaxy is set to the maximum 
possible. The dotted line is the best-fit to the open circles. Right: Probability distribution for Hq. 



This indicates that the possible range in Ho if a reasonable 
B-band galaxy flux is determined for each galaxy will be be- 
tween 21 and 44 km s~^ Mpc~^. In any case, our conclusions 
still stand that too small a value of Hq is determined by this 
model. 



5 DISCUSSION 

Whilst the time-delay distributions are consistent with re- 
processing in a disc, and fluxes and difference spectra can 
be flt by such a model (provided an estimate for the nuclear 
dust extinction is used), the distances determined using this 
method lead to a significantly smaller value of Ho (at least 
a factor of ~1.6) than determined by other methods - those 
that rely on the distance ladde r e.g., the HST Key P roject 
Ho = 72 ± 8 km s"^ Mpc"^ ('Preedma n et al.l 1200 il l, and 
other direct methods such as gravitational le nsing time de- 
lays and the Sunyaev-Zeldovich effect (e.g. [ Schmidt et al.l 
l2004l : iBonamente et al.l l2006l : ISchechted 120051 ). AU the ob- 
servational evidence currently points towards a value of Ho 
that is around 70 km s~^ Mpc~^. As both direct methods 
and distance ladder methods both point to values of Ho 
much greater than what we find, we interpret our result in 
terms of what this large difference in Ho implies for the 
simple disc reprocessing model we have assumed. We also 
note that there is a large scatter in the distances measured 
around the best-fit model. 

We now briefly recap the main assumptions of our 
model, as this is relevant to our discussion. Firstly, we as- 
sume that the variable ionizing radiation travels outwards 
at the speed of light, so that R = cr, and that there are 
no significant time delays other than the light travel time. 
We also assume that the geometry is of a flat surface, and 
that there is azimuthal symmetry, so that the full annulus 
contributes, giving the surface area at radius R as d(Area) 
= 2nRdR. We assume that there is thermal reprocessing, 
so that the relationship between surface brightness, wave- 
length and temperature is like that of a blackbody. Pinally, 
we correct the slope of the difference spectrum to oc u^^^ 



assuming dust extinction from both our Galaxy and instrin- 
sic to the AGN. Por the purposes of the following discussion, 
we list some scaling relations from our model. The AGN 
flux / cx BOc^T^ cos i/D^ , thus the distance to the object 
D oc ct(B6 cosi/ fy^^ , and Ho oc l/D. Here B is the surface 
brightness, and we assume 6 = 2Tr (azimuthal symmetry). 



5.1 Scatter in distances 

In our Hubble diagram (Fig. 1111) . we note that there is a 
large scatter in the distances around the best-fitting cosmo- 
logical model. An extreme example of this is the two objects, 
NGC 7469 and NGC 5548, that have almost identical red- 
shifts, yet, in the Hubble diagram they are separated by ap- 
proximately 3 magnitudes in the distance modulus. Whilst 
some of the scatter seen is likely due to adopting an inclina- 
tion angle of i = 45° for all objects (as might be expected 
for a type I AGN), this effect is only small - the distance 
goes like (cosi)^ ^^. According to unified AGN schemes (e.g. 
lAntonuccilll993l ). for type I objects where we see the broad- 
emission line region, i < 60° and so 0.707 < (cosi)^'^^ < 1. 
Thus, whilst the inclination does cause scatter, it cannot 
account for the large difference between NGC 7469 and 
NGC 5548. Comparing the lightcurves, AGP and CCPs for 
these two objects (Pigs.[T}, it is seen that NGC 5548 is more 
slowly varying, and hence has a broader ACP, suggesting 
that maybe a difference between the two objects could be 
due to biased time delay measurements. 

To investigate this, and further possible causes of this 
scatter. Pig. [T2] examines correlations between the residuals 
of the fit to the distances and various AGN properties and 
parameters of the fit. Correlations with the mass, luminosity, 
fraction of Eddington luminosity, E{B—V), ACF width, and 
fractional contribution of the galaxy to the minimum flux in 
the V were looked for, but importantly, none were found 
(see Pig. [T2j. In this figure, we highlight NGC 7469 (open 
square) and NGC 5548 (open circle). While the width of the 
ACPs differ largely between NGC 7469 and NGC 5548, there 
is no correlation between the ACP width and the residuals. 
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as would be expected if this was the cause of the scatter. 
Similarly, no correlation with E{B — V) suggests that our 
adopted extinction law does not cause the scatter. 

One possible explanation for the residuals could be due 
to the fact that we fix the galaxy in the B-band to be zero, 
as in our method we have no way of determining this. This 
is clearly an incorrect assumption, and better constraints on 
the galaxy contribution may significantly improve matters. 

Alternatively, reprocessing may not be the only source 
of the variability. For instance. X-ray and optical lightcurves 
for NGC 5548 are seen to be correlated, but the fractional 
variabili ty in the optical is larger than those seen in the 
X-rays (|Uttlev et al.ll2003l ). However, this is the long-term 
optical and X-ray variability, and the short-term variabil- 
ity is seen to be consis tent with reprocessing in this object 
l|Suganuma et al.ll2006l '). 

5.2 Systematic offset in distance 

Having found that Hq is at least a factor of 1.6 too small, we 
discuss possible causes for this systematic offset in distance. 
One suggestion for the difference is that the lags measured 
via cross-correlation are biased to high values, hence giv- 
ing high distances (and low Ho)- As the cross-correlation 
function (CCF) is the convolution of the delay map with 
the auto-correlation function (ACF), changes in the driving 
lightcurve properties (i.e. in the ACF), can lead to a differ- 
ent CCF, and hence a different lag, without changing the 
delay map. If this is the case, one might expect to see a 
correlation between the residuals of the Hubble-diagram fit 
and the full- width half maximum of the ACF. However, as 
discussed above, no such correlation is evident in Fig. 1121 

Another possibility is that the irradiating source does 
not illuminate the entire disc. If this is the case, then the 
solid angle for the reprocessed optical emission at each annu- 
lus in the disc, will be less than 2nR dR cos i and hence the 
distances will be closer than determined and Ho would be 
greater. We briefly consider several possible causes of this. 
Firstly, the irradiating source may be anisotropic, only illu- 
minating specific parts of the disc (hence 9 < 2-k). Secondly, 
the surface of the disc may not be fiat enough. For instance, 
if the surface was bumpy, one can imagine that only the front 
side of the bump would be illuminated and that large bumps 
would actually place some of the disc in shadow. For there 
to be a factor of ~ 2 increase in Ho, there would need to be 
a factor of 4 increase in the actual disc fiux (as oc 
and Ho oc 1/D). Therefore one quarter of the disc would 
need to be illuminated, and the remaining three quarters 
not for the true disc flux to be four times larger than that 
observed. In the l amppost model for disc illumination (e.g., 
ICoUin et al.ll2003l ) this would be hard to achieve as in this 
model the irradiating source is high above the central black 
hole and illuminates it in all directions. However, if the inci- 
dence angle is low, so that irradiating photons can graze the 
disk surface, then small undulations in the disk surface, such 
as might arise from spiral density waves, may be sufficient 
to shadow a significant fraction of the disk surface. Alterna- 
tively, if the illumination of the disc was more localised, for 
instance due to magne tic fiares in the dis c corona just above 
the disc surface fe.g.. [Collin et al]l2003h then, it is easy to 
picture that only a small amount of the disc is irradiated. 

It is important to consider whether additional extinc- 
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Figure 12. The residuals of the Hubble diagram fit are compared 
with various AGN properties and parameters of the fit. The open 
circle indicates NGC 5548, and the open square NGC 7469. Here 
we use the reverberation mapping masses for these objects deter- 
mined bv lPeterson et al. 1 1 I2OO41) . 
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tion can be used to explain an increased true flux. As ex- 
tinction is wavelength-dependent, increasing its magnitude 
would also cause a change in shape of the spectral energy 
distribution. Therefore only if the dust grains in the AGN 
are much larger than assumed (i.e. J? > 5.15 that we used 
for the AGN extinction) could more extinction be allowed. 

Furthermore, a scattering medium above the disc (such 
as a disc wind) could slow the outward propagation of the 
irradiating light, by making it more of a random walk dif- 
fusion process, rather than direct path to the disc surface, 
this could increase the mean effective path length so that 
(r) > R/c, and cause incident photons to be distributed 
over a wider solid angle than their direct line of sight. Such 
a process would broaden the transfer function. If scattering 
is important, then the detailed shape of the disk surface may 
be less important. 

Finally, we also consider whether thermal reprocessing 
from dust can account for this offset in distance. iGaskelll 
l|2006t ) has recently proposed that the short-wavelength tail 
of thermal emission from the dusty torus may account for 
the observed increase of optical continuum time delay with 
wavelength. However, the typical size of the dusty torus is 
thought to b e on scales of up to hundreds of liRht days 
llGlassI I2OO4I : iMinezaki etHI l2004l : ISuganuma et al.l |2004 
[2003). The optical tail of the thermal emission from such 
dust would exhibit similar large delays and broadening, af- 
fecting broad structure in the cross-correlation function, but 
not significantly moving the cross-correlation peak that we 
are using here to measure inter-band delays. 

Observationally, we can better understand AGN discs 
by increasing the cadence of multi-band monitoring to im- 
prove the delay measurements and ideally to measure delay 
distributions rather than just mean delays. The main un- 
certainty in our distances arises from uncertainties in the 
delay measurements. Defining the shape of the delay distri- 
bution, for comparison with predictions, would stringently 
test the blackbody disk model, and could detect or rule out 
additional significant broadening mechanisms. 

5.3 Possible cosmological probe? 

A systematic shift in the distance determination that is ap- 
parent in this analysis is not necessarily a problem when de- 
termining cosmological parameters other than Ho- For Qm 
and r^A , it is the shape of the redshift-distance relationship 
that matters. It is therefore possible that even without un- 
derstanding this offset, this method can be used as a cosmo- 
logical probe provided that the systematic shift is redshift, 
and luminosity, independent. As mentioned earlier, AGN are 
numerous throughout the Universe, out to high redshift. As 
they are intrinsically extremely luminous objects, accurate 
photometry can be obtained out to redshifts far greater than 
supernovae type la can be observed, allowing similar accu- 
racy in flu and Q,a with significantly more objects. More 
luminous AGN will need to be observed at greater redshifts, 
however, and although more luminous quasars are known 
to be less variable, quasar variability increases toward the 
blue part of the spectrum (Vanden Berk et al., 20o4 ). From 
studies of 25,000 quasars in the Sloan Digital Sky Survey, 
IVanden Berk et al.l (|2004l ') also find a positive correlation of 
variability amplitude with redshift. The redshift will move 
the observable region of the spectrum into the UV (in the 



emitted-frame) , which is more variable, and likely has less 
galaxy contamination, however the time delays there will 
be shorter. Counteracting this is time-dilation effects, which 
mean the observed time delay will be longer than the emit- 
ted time delay. 

To explore possible constraints on Qm and Qa through 
monitoring programs, we simulate the distances to a sam- 
ple of AGN. With the current dataset we have around 25% 
rms errors in the distance, but the monitoring was only once 
every 2-3 days, improved temporal sampling should signif- 
icantly improve the accuracy of the time-delay measured, 
and hence the distance determined. We therefore assume 
that the distances can be measured to 15% accuracy, which 
is plausible given high-cadence monitoring with a suite of 
robotic telescopes across the globe. To approximate the red- 
shift distribution of variable AGN that might be available 
in practice, we adopted a sample including the 14 nearby 
AGN in this work, augmented by 30 AGN that were discov- 
ered in the magnitude-limited QUESTl v ariability survey 
(Qua sar Equatorial Survey Team, Phase 1, iRengstorf et al.l 
'2004'). This provides known variable AGN over a range of 
redshifts. In the simulation, we determine the distances from 
the redshift, assuming Hm = 0.3, flA ~ 0.7 and Ho = 72 km 
s~^ Mpc~^. We add a random Gaussian-distributed error to 
each of the points, with a mean of 0.0, and 1-a = 15%. We 
then fit Eq[TS]by minimizing x^- 

The results of this simulation are shown in Fig. 1131 For 
comparison with our results, recent constraints from super- 
nov ae observations can be seen for example in Fig. 12 of 
Wo od-Vasev et al] (|2007f ) . For each value of f2 m and Qa op- 
timal scaling is used to determine the best value of Ho- The 
best-fitting model for a fiat Universe {fltot ~ ^Im+^a ~ 1.0) 
is indicated in each figure. The probability distribution for 
Ho in these figures is determined for a flat Universe and also 
when there is no constraint on Qtot- In the case of assuming 
a fiat Universe, for each value of Ho the best-fitting model 
and the probability are determined - the only free parame- 
ter is Qm, as the condition that Qtot = 1.0 determines Qa- 
Clearly, the input model is retrieved from the simulations, 
showing that this method has potential as a cosmological 
probe, if the systematics can be understood, or are redshift 
independent. 



6 CONCLUSIONS 

We have fitted the wavelength-dependent time delays and 
optical SEDs of 14 AGN using a disc reprocessing model, 
which has allowed a measure of the nuclear reddening in 
these AGN, as well as a measurement of the distances to the 
AGN. However, the distances calculated using this method 
imply Ho = 44 ± 5 km s~^ Mpc~^, a factor of 1.6 less 
than the value that all other methods seem to be pointing 
to. We have discussed the basic assumptions of the black- 
body disc model and a number of possible systematic effets 
to which the simple blackbody disc model may be vulner- 
able. Even with a systematic shift, it may be possible to 
use this method to probe Qm and Qa, as it is the shape of 
the redshift-distance relationship that matters. We have pre- 
sented a simulation showing possible constraints from long- 
term monitoring of 44 AGN. 
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Figure 13. Simulation of constraints on Hq, Qm and Qa from 44 AGN. Top ie/t; Distance modulus vs. redshift for the 44 AGN. 
Various cosmological models are shown, with the solid line indicating the best-fitting flat cosmology. Bottom left: Magnitude difference 
between the distance modulus and an empty Universe. Top right: Probability distribution for Hq. The solid line indicates the probability 
distribution when a flat cosmology is assumed, and the dashed line shows the distribution with no constraint on Cltot-Bottom right: 
Probability distribution for ^Im and Ha- Contours indicate 1, 2 and 3-cr confidence limits. 
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